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The Turaev-Viro invariants are scalar topological invariants of compact, orientable 3-manifolds. 
We give a quantum algorithm for additively approximating Turaev-Viro invariants of a manifold 
presented by a Heegaard splitting. The algorithm is motivated by the relationship between topolog- 
ical quantum computers and (2 -I- 1)-D topological quantum field theories. Its accuracy is shown to 
be nontrivial, as the same algorithm, after efficient classical preprocessing, can solve any problem 
efficiently decidable by a quantum computer. Thus approximating certain Turaev-Viro invariants of 
manifolds presented by Heegaard splittings is a universal problem for quantum computation. This 
establishes a novel relation between the task of distinguishing non-homeomorphic 3-manifolds and 
the power of a general quantum computer. 
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The topological quantum computer is among the most 
striking examples of known relationships between topol- 
ogy and physics. In such a computer, quantum infor- 
mation is encoded in a quantum medium on a 2-D sur- 
face, whose topology determines the ground space de- 
generacy. Surface deformations implement encoded op- 
erations. Topological quantum computers are universal, 
i.e., can implement arbitrary quantum circuits. It is nat- 
ural to try to identify the topological origin of this com- 
putational power. 

One answer is that the power stems from the un- 
derlying (2 -I- 1)-D topological quantum field theory 
(TQFT) [1]. The TQFT assigns a Hilbert space to a 
2-D surface S, and a unitary map U{f) : H-^ — >■ Hj:' to 
every diffeomorphism / : S ^> S', subject to a number 
of axioms [2]. However, this answer is not fully satis- 
factory; the definition of a TQFT is involved, and uses 
mathematics that appears in similar form in the theory 
of quantum computation. A second answer, arising in [3- 
6], is that quantum computers' power comes from their 
ability to approximate the evaluation, at certain points, 
of the Jones polynomial of the plat closure of a braid. 

Here we give an alternative topological description 
of the power of quantum computers, in terms of the 
Turaev-Viro 3-manifold invariants. Observe that restrict- 
ing TQFTs to closed manifolds results in scalar invari- 
ants. We show that approximating certain such invari- 
ants is equivalent to performing general quantum compu- 
tations. That is, we give an efficient quantum algorithm 
for additively approximating Turaev-Viro invariants, and 
conversely we show that for any problem decidable in 
bounded-error, quantum polynomial time (BQP), there 
is an efficient classical reduction to the Turaev-Viro in- 
variant approximation problem. The classical procedure 
outputs the description of a 3-manifold whose certain 
Turaev-Viro invariant is either large or small depending 
on whether the original BQP algorithm outputs 1 or 0. 



Turaev and Viro [7] defined a family of invariants for 
compact, orientable 3-manifolds. The original defini- 
tion parameterized the invariants by the quantum groups 
SU(2)fc, for k G N, but it was extended by Barrett and 
Westbury [8] to give an invariant for any spherical ten- 
sor category C. Any compact 3-manifold M is homeo- 
morphic to a finite collection of tetrahedra glued along 
their faces [9]. Beginning with such a triangulation, as- 
sign a certain rank-six tensor F to each tetrahedron and 
a certain gluing tensor d to every edge. The invariant 
TVc(M) is the contraction of the tensor network, which 
can be written out as 

TVc(M)=p-W J2 I[de U (1) 

labclings edges tetrahedra V m n 

if C is multiplicity-free. Here, the sum is over edge la- 
belings of the triangulation by particles from the cate- 
gory C. The index e on d is the label of an edge, while 
the indices i, . . . ,n are the labels of the six edges involved 
in a tetrahedron, ordered and oriented following certain 
rules. The fusion tensor F, the quantum dimensions d 
and the total quantum dimension V are parameters of C. 
\V\ is the number of vertices of the triangulation. The 
topological invariance of TVc (M) follows from the fact 
that any two triangulations of M can be related by a fi- 
nite sequence of local Pachner moves [10], under which 
the above quantity is invariant. In this paper we con- 
sider multiplicity-free unitary modular tensor categories, 
which include the SU(2)fc case, but are not as general as 
spherical tensor categories. 

To formulate a BQP-complete problem [11] of estimat- 
ing the Turaev-Viro invariant, we require a presentation 
of 3-manifolds known as a Heegaard splitting. Consider 
two genus-g handlebodies (e.g., the solid torus for g — 1). 
They can be glued together, to give a 3-manifold, us- 
ing a self-homeomorphism of the genus-^ surface. The 
set of orientation-preserving self-homeomorphisms mod- 
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FIG. 1: A Dehn twist is a 2n rotation about a closed curve. 
The Dehn twists about the Sg—1 curves shown above generate 
the full mapping class group of the genus-gr surface [12]. 

ulo those isotopic to the identity form the mapping class 
group MCG(g) of the surface. It is an infinite group gen- 
erated by the 3g — 1 Dehn twists illustrated in Fig. 1. 
A Heegaard splitting thus consists of a natural num- 
ber g and an element x € MCG{g), defining a manifold 
M{g,x). Every compact, orientable 3-manifold can be 
obtained in this way, up to homeomorphism. 

Theorem 1 For any fixed multiplicity-free unitary mod- 
ular tensor category C , there is a quantum algorithm that, 
given S,e > 0, g £ N and a length-m word x in the 
Dehn-twist generators of MCG{g) from Fig. 1, runs in 
fime poly (g, m, log 1/5, 1/e) and, except with probability 
at most 5, outputs an approximation ofTYc{M{g,x)) to 
within ±X'2(ff-i) 

Conversely, for C the category associated to SU(2)fe or 
S0(3)fc for fc > 3 such that k -\- 2 is prime, it is BQP- 
hard to decide whether 'D'^^^~3^ TYc{M{g, x)) is greater 
than 2/3 or less than 1/3. More precisely, given any 
quantum circuit TofT two-qubit gates acting on n qubits 
|0"), with output either or 1, one can classically find 
in polynomial time a word the standard 

Dehn-twist generators of yiCG{g), with g = n -\- 1 and 
m = poly(r), such that 

|Pr[T outputs 1] -^^(i-ff) TYc{M{g,x))\ < 1/6 . (2) 

The additive approximation error is exponential in g. 
Complexity-theoretic reasons make it unlikely that quan- 
tum computers can efficiently obtain a multiplicative or 
otherwise presentation-independent error [13]. 

In fact, a similar statement to Theorem 1 also holds for 
approximating the Witten-Reshetikhin-Turaev (WRT) 
invariants [14, 15]. For any g, a modular category C 
can be used to define a projective representation pc.g '■ 
MCG{g) — GL{'Hc.g)- This representation will be given 
below. The WRT invariant for a 3-manifold M{g,x) is 
then given by a matrix element 

WRTciM{g,x))=Va-' {vcj pc,g{x) \vc,g} , (3) 

where \vc.g) G 'Hc,g is a certain unit-normalized vector. 
As the representation is projective, WRTc is a 3-manifold 
invariant only up to a multiple of e,^'^'^'^/^^ where c is 
called the central charge. (Eq. (3) is the Crane-Kohno- 
Kontsevich presentation [16-18] of the WRT invariant, 
which is more commonly defined in terms of a Dehn 



surgery presentation of M. Equivalence of these defi- 
nitions for C = SU(2)fc is shown in [19]; see also [20, 
Sec. 2.4].) 

The fact that Eq. (3) indeed gives an invariant can be 
established by studying the problem of when two Hee- 
gaard splittings and {g',x') describe homeomor- 
phic manifolds. Since taking the connected sum of a 
manifold M with the 3-sphere does not change the 
manifold, i.e., M^S^ = M, the standard Heegaard split- 
ting of into two genus-one handlebodies allows defin- 
ing a "stabilization" map {g,x) i— (5 + 1,5) such that 
M{g, x) = M{g -\- 1, x). A general theorem of Reidemeis- 
ter [21] and Singer [22] asserts that M{g,x) = M{g',x') 
if and only if {g,x) and {g',x') are equivalent under sta- 
bilization and the following algebraic equivalence relation 
for the case of equal genus [23] 

{g,x) = {g,x') if a; = with y,z e MCG+(g) . (4) 

Here MCG+(5) C MCG(g) is the subgroup of self- 
homeomorphisms (classes) of the genus-g surface that 
extend to the genus-g handlebody. Invariance of 
WRTc(M(g,x)) now follows essentially from the fact 
that \vc,g) is invariant under the action of MCG'^{g). 
The Turaev-Viro and WRT invariants are related by 

TVc(Af) = [WRTc (M) I' (5) 

as shown by Turaev [24] and Walker [2] (see also [25, 26]). 
In [27], Eq. (5) is discussed in the category-theoretic 
formalism used here. Identities (3) and (5), together 
with density and locality properties of the representa- 
tions pc,g, are the basis of our BQP-completeness proof. 

Previously, a quantum algorithm for approximating 
the SU(2)fc Turaev-Viro and WRT invariants was given 
by Garnerone et al. [28], assuming the manifold is spec- 
ified by Dehn surgery rather than a Heegaard splitting. 
BQP-hardness of the approximation was left as an open 
problem. In unpublished work, Bravyi and Kitaev have 
proven the BQP-completeness of the problem of approx- 
imating the SU(2)4 WRT invariant of 3-manifolds with 
boundary [29], where the manifold is specified using 
Morse functions. We remark that one can use Arad and 
Landau's quantum algorithm for approximating tensor 
network contractions to compute the Turaev-Viro invari- 
ant of a triangulated manifold [30] . While this algorithm 
would run polynomially in the number of tetrahedra, its 
precision depends on the order in which tensors are con- 
tracted and may be trivial. 

We will only briefly describe the space "Hc.g, the rep- 
resentation pc,g : MCG(g) — >■ GL{'Hc,g) and the state 
\vc,g) G ^C,s from Eq. (3). Details are in [16-18, 27]. 

Let C be a multiplicity-free unitary modular tensor 
category. It specifies a set of particles i with quantum 
dimensions di > 0, and including a trivial particle 0. 
The total quantum dimension is V — "^l- ^ addi- 

tionally specifies a particle duality map i t-^ i* , fusion 
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FIG. 2: Three examples of decompositions of the genus-two 
surface E2 into three-punctured spheres. In each case, a triva- 
lent adjacency graph of the punctured spheres is shown in red. 



unitary, with one control qudit, and its action is given by 



(7) 



Most presentations of modular tensor categories do not 
explicitly provide values for S'jj,. However, as discussed 
in [2], S'k can be calculated by the identity 




rules, F-symbols Flj™ and i?-symbols Rj''. These ten- 
sors obey certain identities, such as the pentagon and 
hexagon equations, which can be found in, e.g., [27, 31]. 

Let g G N, g > 2. The space Hcg can be defined by 
specifying an orthonormal basis. Decompose the genus-5 
surface Sg into three-punctured spheres (or "pants") by 
cutting along 85 — 3 noncontractible curves, as illustrated 
in Fig. 2. Dual to such a decomposition is a trivalent 
graph r. Direct arbitrarily the edges of T. A basis vec- 
tor |€)p is a fusion-consistent labeling of the edges of F by 
particles of the category C. Fusion-consistency is defined 
by the fusion rules, i.e., a set of triples (z, j, k) that are al- 
lowed to meet at every vertex, and particle duality, which 
switches the direction of an edge, replacing a label i by 
the antiparticle i* . Define the states Br {|^)r}^ to be 
orthonormal, and their span to be Hcg- Note that this 
definition gives a natural encoding of Hcg into qudits, 
with one qudit to store the label of each edge of F. The 
directed graph F can be stored in a classical register. 

The above definition depends on F, but alternative 
pants decompositions simply represent different bases Br 
for the same Hilbert space. To convert between all possi- 
ble pants decompositions of Eg we need two moves, each 
corresponding to a local unitary operator. 

The F move relates bases that differ by a "fiip" of a cut 
between two three-punctured spheres. In the qudit en- 
coding, it is a five-qudit unitary, with four control qudits. 
Its action is given by 




(6) 
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The S move applies when two boundaries of a single 
three-punctured sphere are connected. It is a two-qudit 



(8) 

(The last expression uses ribbon graph notation.) 

The action pc^g of MCG(5) on Hcg can now be speci- 
fied by the action of the Dehn-twist generators on basis 
vectors. For a Dehn twist about a curve ct, apply a se- 
quence of F and 5" moves to change into a basis By, i.e., a 
pants decomposition of Eg , in which a divides two three- 
punctured spheres. In such a basis, the Dehn twist acts 
diagonally: if the edge of F crossing a has label i, the 
twist applies a phase shift of . 

To complete the definition of WRTc(M((7, a;)) from 
Eq. (3), it remains to define the state \vc,g)- As on the 
right-hand side of Eq. (7), decompose Eg with a merid- 
ional cut through each handle. Then |wc,g) is the state in 
which every edge of F is labeled by 0, the trivial particle. 

Let us now prove Theorem 1. Although not obvious 
from Eq. (1), the original tensor-network-contraction- 
based definition of the Turaev-Viro invariant. Theorem 1 
is a straightforward consequence of the definition based 
on the representation pc,g, and of known density results. 

The Turaev-Viro and WRT invariants for M{g,x) can 
be approximated essentially by implementing pc^g{x). 
The algorithm maintains a classical register storing the 
graph F, together with a quantum register containing the 
current state in Hcg in the basis Br- If C has N parti- 
cle types, the algorithm uses an A'^-dimensional qudit for 
each edge of F. Then pc^g{xj) can be applied by using a 
sequence of F and S moves, i.e., certain local unitaries, 
to change to a basis in which Xj acts diagonally. Since 
Xj is one of the generators from Fig. 1, starting with the 
graph F of Fig. 3(a) (for which every edge is labeled 
in \vc.g)) at most one F and one S move are needed. 
An estimate to within e of the desired matrix element 
{vc,g\ pc.g{x)\vc,g) can be given, except with probabil- 
ity 5, using 0(log(l/5)/e^) Hadamard tests, as in [3]. 

To prove BQP-hardness we reduce from the BQP- 
complete problem of deciding whether |(0^|T|0^)p is 
larger than 5/6 or less than 1/6, given the 5-qubit quan- 
tum circuit T [3]. Let C be the modular tensor category 
associated with SU(2)fc or S0(3)fc, with A: > 3 and k + 2 
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FIG. 3: (a) A p-qubit state \z), z G {0, 1}*, can be encoded 
into Hc,g for the genus-g handlebody. (b) Any two-qubit gate 
can be approximated within the codespace using the Dehn 
twists involving the two corresponding handles. 



prime. Given T consisting of T two-qubit gates, our aim 
is to construct efficiently the Heegaard splitting (g, x) of 
a manifold M = M{g, x) such that P^Ci-s) TVc(M) ap- 
proximates 1(0^1 T |0^)p. As illustrated in Fig. 3(a), we 
use one handle of a genus-g handlebody to encode each 
qubit. Such a labeling is fusion-consistent, and the en- 
coding of the initial state |0^) is exactly \vc,g) G '^c,g- 
As shown in [32, 33], for C = 80(3)^ the representa- 
tion pc^g has a dense image, up to phases, in the group 
of unitaries on Hcg, for g > 2. By the density for 
g = 2 and the Solovay-Kitaev theorem [34] , therefore any 
two-qubit gate can be approximated in the codespace to 
precision 1/(6T) by applying a (log r)'^(^)-long sequence 
of the five Dehn twists shown in Fig. 3(b). This holds 
also for C — 811(2)^, as 80(3)^ is just the restriction 
of 8U(2)fe to particles with integer spins. Thus we ob- 
tain a polynomial-length word x = xi - ■ ■ Xpoiy(T) in the 
Dehn-twist generators whose action approximates T on 
the codespace. Then {vc,g\ pc,g{M{g,x))\vc,g) approxi- 
mates (O^l T [Qf). 

This work demonstrates how quantum physics, in the 
form of TQFTs, can inspire new quantum algorithms for 
problems based on topology and tensor networks. The 
approach taken here realizes in a sense the traditional 
vision of quantum computers as universal simulators for 
physical systems, but with a different outcome: it pro- 
vides a purely mathematical problem whose difficulty ex- 
actly captures the power of a quantum computer. 
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